Half-metallicity in NiMnSb: 
a Variational Cluster Approach with ab-initio parameters 
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Electron correlation effects in the half-metallic ferromagnet NiMnSb are investigated within a com- 
bined density functional and many-body approach. Starting from a realistic multi-orbital Hubbard- 
model including Mn and Ni-d orbitals, the many-body problem is addressed via the Variational 
Cluster Approach. The density of states obtained in the calculation shows a strong spectral weight 
transfer towards the Fermi level in the occupied conducting majority spin channel with respect to 
the uncorrelated case, as well as states with vanishing quasiparticle weight in the minority spin gap. 
Although the two features produce competing effects, the overall outcome is a strong reduction of 
the spin polarisation at the Fermi level with respect to the uncorrelated case. This result emphasizes 
the importance of correlation in this material. 



I. INTRODUCTION 

More than twenty years ago de Groot et al} carried 
out electronic structure calculations for the half-Heusler 
compound NiMnSb which showed peculiar magnetic fea- 
tures leading to the discovery of a new class of materials, 
the so-called half-metallic ferromagnets. Such materials 
differ from conventional ferromagnets in that they dis- 
play a gap in one of the spin channels only. The con- 
cept of half metallicity boosted the research in spintron- 
ics - an emergent technology which makes use of spin 
and charge of electrons at the same time. Spintronic 
applications such as spin- valves, polarized electron in- 
jectors/detectors or devices using tunneling and giant 
magneto-resistance effects promise to revolutionize mi- 
croelectronics once highly polarized electrons can be in- 
jected efficiently at room temperatures^i^,. 

Unfortunately, the theoretically predicted ideal full 
spin polarization of half-metals has not yet been found 
experimentally. As a matter of facts, experiments show 
that full polarization is lost at temperatures of the order 
of room temperature, and even at lower temperatures dif- 
ferent factors such as structural inhomogeneities, as well 
as surfaces and interface properties may suppress \\^. 

Despite the fact that high quality NiMnSb films have 
been successfully grown, they were not found to repro- 
duce the half-metallic character of the bulk suggested 
by spin-polarized positron-annihilation^. Values of spin 
polarization were reported between 40% in spin-resolved 
photoemission measurements® up to 58 ± 2.3% by super- 
conducting point contact measurements at low temper- 
atures^ (see also Refs. lolfTol). The discrepancy between 
theoretical calculations-'^ and the above mentioned exper- 
imental facts were attributed to surface and interface ef- 
fects. Consequently, different surface and interfaces of 
NiMnSb were theoretically investigated by de Wijs and 
de Groot^ , which demonstrated that half-metallicity can 
be preserved at the surface and/or interface by suitable 
reconstruction^. The theoretical situation is complicated 



by the fact that the spin polarisation (or, more precisely 
the tunneling magnetoresistance) displays a substantial 
uniaxial anisotropy in this materialj^ 

Recently, finite-temperature correlation effects were 
addressed in several half-metals^ ' -'^^'^'^d^'^^d^ . For 
NiMnSb, a Local Density Approximation plus Dynam- 
ical Mean Field Theory calculation (LDA-f DMFT>i^ 
showed the appearance of so-called non-quasiparticle 
(NQP) states. These states originate from spin-polaron 
processe s ^^d® , whereby the spin-down low-energy elec- 
tron excitations, which are forbidden for half- metallic 
ferromagnets in the one-particle picture, turn out to be 
possible as superpositions of spin-up electron excitations 
and virtual magnon o^d'^d^ . Here, we extend this study 
by adopting the Variational Cluster Approach (VGA), 
which includes correlations beyond the locality captured 
by DMFT. In addition, the VGA is based on exact di- 
agonalisation, which is more appropriate than the dia- 
grammatic method (FLEX) adopted in Ref. [Ij to solve 
the impurity problem. In a previous paper we used 
the VGA to investigate the spin polarization in NiMnSb 
taking into account only the Mn-d orbital basis set. Our 
calculations showed that the Mn-d-only basis set is not 
sufficient to appropriately describe the low energy spec- 
trum of NiMnSb around the Fermi level. For this reason, 
in the present work we adopt a multi-orbital Hubbard- 
typc Hamiltonian which includes all 10 Mn and Ni-d or- 
bitals. Our present calculation confirms that the inclu- 
sion of the latter is essential for a proper description of 
ferromagnetic properties and of the minority spin gap in 
NiMnSb. 

Our results support the existence of states within 
the minority spin gap in agreement with previous 
LDA-I-DMFT calculations'^-^^. In addition, they indi- 
cate that these so-called nonquasiparticle states indeed 
have a vanishing quasiparticle weight at the Fermi en- 
ergy. At the same time, our results predict a correlation- 
induced spectral weight transfer for the majority spin 
states. The combination of these two effects yields a 
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polarization whose energy dependence is in qualitative 
agreement with experiments. These calculations lead to 
the conclusion that even in the presence of medium-size 
interactions, electron correlations significantly affect the 
spin polarisation in half- metals. 

This paper is organized as follows: in section |TT] we 
present the methods used to investigate the electronic 
structure of NiMnSb. In particular, in Sec. Ill Al we 
describe the ab-initio construction of the many-body 
model Hamiltonian. Specifically, the uncorrelated part 
of the Hamiltonian for excitations in the vicinity of the 
Fermi level is obtained from the so-called downfolding 
technique2£i2i within the Nth-order muffin tin orbital 
(NMTO) method. In Sec. IIIBl we give a short sum- 
mary of the VGA approach. We present and discuss our 
results in Sec. lIIII In particular, in Sec. lIII Al we evaluate 
the density of states within VGA and discuss the results 
in the framework of previous calculations. In Sec. IIII Bl 
we discuss fc-dependent spectral properties, namely the 
spectral function and the self-energy. Finally, spin polar- 
ization and its comparison with experiments is discussed 
in Sec. lIII Cl and the summary of the results is presented 
in Sec. |TVl 



II. ELECTRONIC STRUCTURE 
CALCULATIONS FOR NIMNSB 

The intermetallic compound NiMnSb crystallizes in the 
cubic structure of MgAgAs type (Clb) with the fee Bra- 
vais lattice (space group FA3m = Tf). This structure 
can be described as three interpenetrating /cc lattices of 
Ni, Mn and Sh with the lattice parameter a — 11.20ao 
(oq = Bohr radius), respectively. The Ni and Sb sublat- 
tices are shifted relative to the Mn sublattice by a quar- 
ter of the [111] diagonal in opposite directions, see also 
Fig. [TJ The important aspectaii^i^^i^^i^iS^ which deter- 
mine the behavior of electrons near the Fermi level, as 
well as the half-metallic properties are the interplay be- 
tween the crystal structure, the valence electron count, 
the covalent bonding, and the large exchange splitting of 
Mn-d electrons. For the minority spin gap opening, not 
only the Mn-d-Sb-p interactions, but also Mn-d-Ni-d in- 
teractions have to be taken into account. In addition, the 
loss of inversion symmetry produced by the Clh structure 
(the symmetry lowering from Oh in the L2i structure to 
Td in the Clb structure) are important for these effects. 
The existence of sp-valent Sb is crucial to provide stabil- 
ity to this compound. 

The crystal structure is shown in Fig. [TJ The positions 
occupied by atoms are represented by spheres. For illus- 
trative purposes, in Fig. [1] the radii of the spheres were 
chosen arbitrarily. The actual muffin-tin radii used in 
the calculations are R^lj, = 2.584; = 2.840; Rflj. = 
2.981 and Rfij' — 2.583 (atomic units) for the vacant po- 
sition situated in (1/4, 1/4, 1/4). The LMTO-ASA basis 
used for the self-consistent calculations contains the spd- 
partial waves for Mn and Ni, the spfrf/j-partial waves 



for Sb and sfpd^-partial waves for the empty sphere E. 
(I) means that the Z-partial waves are downfolded within 
the selfconsistent calculations. 




FIG. 1: (Color online) The conventional unit cell for the semi- 
Heusler NiMnSb compound: Sh (large, blue spheres) and Mn 
(medium, green spheres) sit on the same faces of the large 
cube (shown with thin, solid lines), Ni (small, yellow spheres) 
forms a separate (small) cube drawn using dashed lines. In 
addition, the atoms belonging to the reference system of the 
VGA calculation are connected by a thick solid line. Sb or- 
bitals are downfolded and, therefore, not included directly in 
the model calculation. 



A. Ab-initio construction of the model 
Hamiltonian 

In order to construct the effective low-energy Hamilto- 
nian to use in our VGA calculation, we employed the Nth 
order mufhn-tin-orbitals scheme within the LMTO-ASA 
basis set26. rpj^^ NMTO method^^'-Si can be used to gen- 
erate truly minimal basis sets with a massive downfolding 
technique. Downfolding produces minimal bands which 
follow exactly the bands obtained with the full basis set. 
The minimal set of symmetrically orthonormalized NM- 
TOs is a set of Wannier functions. In the construction of 
the NMTO basis set the active channels are forced to be 
localized onto the eigenchannel RZm, making the NMTO 
basis set strongly localized. 

Fourier-transformation of the orthonormalized NMTO 
Hamiltonian, H^^^{\i), yields on-site energies and hop- 
ping integrals, 

in a Wannier representation, where the NMTO Wannier 
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functions 



^Rm 



are orthonormal. Here, t 



A-B,R'-Il 



de- 



notes the hopping term from orbital m of atom B on site 
R to the orbital m' of atom A on site R' {A and B are ei- 
ther Ni or Mn) Further information concerningtechnical 
details of the calculation can be found in Ref.T5ll27l. 

In a previous paper^*" we discussed the chemical bond- 
ing and computed model Hamiltonian parameters for the 
semi-heusler NiMnSb using only Mn-d Wannier orbitals. 
As mentioned above, not only the Mn-d-Sb-p, but also 
Mn-d-Ni-d interactions are required to open a gap in 
the minoriy spin channel: the minority occupied bond- 
ing states are mainly of Ni-d character, while the unoc- 
cupied anti-bonding states are mainly of Mn-d character. 
Therefore, in the present work we consider an enlarged 
NMTO-basis consisting of Ni- and Mn-d orbitals which 
span an energy window of about ±3eV around the Fermi 
energy. 



The matrix elements for the on-site energies 



A-A,0 



are 



m,m — given 
which m — 1, • • ■ , 



by (we use the 
, 5 corresponds to 
-y^} in the order) 



convention in 
the d orbitals 



(-1411, -1411, -1411, -721, -721) , (2) 



ef^ = (-2439, -2439, -2439, -2679, -2679) 



(3) 



The nearest-neighbour hopping terms are given by (Ai 
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and the next-nearest-neighbour terms (A2 = (i, — i, 0)) 
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Here, all hoppings are given in units of meV, and only one 
representative hopping integral is shown for each class. 
Other hopping terms can be derived from proper uni- 
tary transformation using crystal symmetry (see, e.g., 
Ref. m for details). As one can see, the largest hoppings 
occur between the Wannier orbitals located on Ni- and 
Mn atoms. In addition, there are further hopping terms 
in the Hamiltonian, which we don't show here for sim- 
plicity. We have taken into account hoppings up to a 



range of r = 2.0a. Neglected hoppings are about a fac- 
tor 30 smaller than the largest nearest-neighbor hopping. 
The non-interacting part of the effective Hamiltonian for 
NiMnSb, thus, has the form 



R',R,cr {A,B,7n',m} 



,A-B,R'-R t 



(6) 



To take into account correlation effects, we add the spin- 
rotation invariant interaction Hj 



H^-\ E E 

R,A,cr,CT' m,n,o,p 



(7) 

In Eqs. and ([7]), CARma {cARma) ^I'e the usual 
fermionic annihilation (creation) operators acting on an 
electron with spin a at site R in the orbital m of atom A. 
The Hamiltonian we are using includes spin- and pair-flip 
terms, as especially spin flip processes are important for 
a correct description of non-quasiparticle states^ii^. For 
a realistic description of Coulomb interactions, the ma- 
trix elements Umnop can be computed for the particular 
material in terms of effective Slater integrals and Racah 
or Kanamori coefhcients^i2S. We used for both Mn and 
Ni the following effective Slater parameters: ^^M^yNj = 
1.26eV, F^^^^^. = 5.58eV and F^^^^. = 3.49eV, which 
correspond to U = Ummmm = 2.0eV and give an average 
J — 0.65eV. In addition, we checked that our results do 
not depend significantly on the chosen U and J-values 
by performing additional calculations for U — 2.5eV and 
U = 3eV with J — 0.77eV. The range of values corre- 
sponds to the one used in previous work s^^i^^'^^i'^^ . 

The on-site energies calculated in NMTO already con- 
tain effects from the Coulomb interaction at the LDA 
mean-field level. While this double counting can be ab- 
sorbed into the chemical potential when only one set 
of degenerate orbitals is taken into account (see e.g. 
Refs. mIssI ) this is generally not possible if the full d- 
orbitals and/or the orbitals of more than one inequiva- 
lent atom are used as basis set. Consequently, in order to 
avoid a double counting of the Coulomb contribution, the 
corresponding Hartree terms have to be subtracted from 
HqM- This is achieved by replacing the onsite energies 
e± obtained from NMTO witb^ 



(n^) denotes the occupation of the NMTO band asso- 
ciated with atom A and orbital m. Notice that this is 
different from the occupation of the orbital m in atom 
A. The present double-counting procedure corresponds 
to the "around-mean-field" scheme^ii^. In order to test 
the dependence on the double-counting schemes, we per- 
formed calculations considering different schemes, such as 
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the "fully-localized" scheme and a combination of these 
two schemes as discussed in Ref.lSJ. From these tests, we 
find that our LDA+VCA results for the spectral function 
remain unchanged for energies within Ep ±leY, while for 
energies outside this range (from ±leV to about ±3eV) 
a redistribution of Mn- and Ni- states is obtained. 



B. Variational Cluster Approach 

To solve the many-body Hamiltonian ©-I- ([7]), we 
employ the Variational Cluster Approach ■^^i^. This 
method is an extension of Cluster Perturbation Theory 
(CPT)^^>ii. In CPT, the original lattice is divided into 
a set of disconnected clusters, and the inter-cluster hop- 
ping terms are treated perturbatively. VCA addition- 
ally includes "virtual" single-particle terms to the cluster 
Hamiltonian, yielding a so-called reference system, and 
then subtracts these terms perturbatively. The "opti- 
mal" value for these variational parameters is determined 
in the framework of the Self-energy Functional Approach 
(SFA)^i^, by requiring that the SFA grand-canonical 
potential f2 be stationary within this set of variational 
parameters. In this work, we only include the chemical 
potential of the cluster as a variational parameter, which 
is necessary in order to obtain a thermodynamically con- 
sistent particle density-- . It is not necessary to include 
a ferromagnetic field in order to obtain a ferromagnetic 
phase, since the symmetry can be broken already at the 
finite-cluster level^. In this paper, we use a new method, 
described in Ref.H^, to carry out the sum over Matsubara 
frequencies required in the evaluation of whereby an 
integral over a contour lying a finite (not small) distance 
from the real axis is carried out. 

As a reference system we use a cluster of two sites, 
representing one Mn-atom and one Ai-atom (see Fig.[T]), 
each having the fuU-d manifold of five orbitals. Since 
we have to consider all five orbitals for each atom, it 
is very difficult to use larger clusters, which have to be 
exactly diagonalized many times in combination with the 
variational procedure. 



III. RESULTS 

A. Density of states 

In order to study the influence of correlations on the 
half-metallic gap, we first display the spin-resolved local 
density of states in Fig.[2l Here, we present a comparison 
of the results obtained from LSDA with the results from 
our LDA-I-VCA calculation. The LSDA-DOS is mainly 
characterized by a large exchange splitting (about 3eV) 
of the Mn-d states, leading to large spin moments on the 
Mn-site (3.72/xb). A small induced ferromagnetic mo- 
ment is present on Ni (0.29/zb), while the Sb moment 
(O.OG/is) is anti-parallel to the Mn moment. Overall, the 



calculated moments are in very good agreement with pre- 
vious ab-initio resultsii^i^iS^i^l. The existence of large 
localized Mn moments of about 3.78/xb has been verified 
experimentally by neutron diffraction''^ as well as by the 
sum rule of the x-ray magnetic circular dichroism spec- 
tra*^ . These two experiments also confirm the magnitude 
of the LSDA-computed moments for Ni and Sb. The gap 
in the minority spin channel is about 0.5eV wide and 
the total magnetic moment has an integer value of 4//^. 
Note, that in Fig.[2]in LSDA the minority occupied bond- 
ing states are mainly of Ni-d character, while unoccupied 
anti-bonding states are mainly of Mn-d character. It was 
pointed out' that the opening of a gap is assisted by 
5*6 through the symmetry lowering with the consequence 
that the distinction between Mn-t2g and Sb-p character 
of the electrons is lost. 

Concerning the LDA-I-VCA results, we find a total 
magnetic moment of 3.7/iB, which is in reasonable agree- 
ment with experimental values^"^. In direct comparison 
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FIG. 2: (Color online) Density of states for NiMnSb obtained 
from LDA-fVCA (red, solid line) for values of the average 
Coulomb and exchange parameters U — 2eV and J — 0.65eV 
for both Mn and Ni atoms, in comparison to results obtained 
from LSDA (black, dashed line). 

to LSDA, our results show that correlation effects do not 
affect too strongly the general picture of the minority spin 
DOS for energies which are more than 0.5 eV away from 
the Fermi energy. In the range 0.5eV ^ E — Ep < 3eV, 
unoccupied Mn- states are visible in the minority spin 
sector, similarly to LSDA. However, these states are 
shifted to lower energies due to a slight reduction of the 
Mn- exchange splitting generated by the many-body cor- 
relations. Just above the Fermi level, NQP states are 
present, with a peak around the energy of 0.06eV. It is 
important to note that these states were also obtained in 
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previous calculations using a LDA+DMFT many-body 
approacbi^ at finite temperatures. In comparison to the 
DMFT description, the non-local correlations captured 
by VGA enhance the spectral weight of the NQP states 
and slightly shift their position. This fact leads to the 
conclusion that even a local DMFT description is suf- 
ficient to demonstrate the existence of NQP states as 
discussed previously The spectral weight of the NQP 
states is large enough, so that we expect them to be well 
pronounced in corresponding experimental data. While 
model calculations for single-band Hamiltonians^ suggest 
that NQP states should only touch the Fermi level with 
zero weight at T = OK, in our VGA calculation they 
maintain a finite weight at the Fermi level, thus leading 
to a reduction of spin polarisation, even at T = OK. In 
the LSD A- results, the bonding states below the Fermi 
level have dominant Ni-d character and are responsible 
for the gap formation. While these states form a single 
peak at -1.5eV in LSDA, the LDA+ VGA-results show 
a splitting into two peaks centered around — leV. One 
of these peaks is pushed closer to the Fermi level, while 
the other one is shifted to higher energies. The latter 
correlation effect is also seen in previous LSDA-I-DMFT 
results^^. 

A significantly stronger effect caused by many-body 
correlations is visible in the majority-spin channel (see 
Fig. (2) . Here we discuss the behaviour in the same en- 
ergy range within ±3eV around the Fermi level, since 
this is the energy window spanned by our NMTO ba- 
sis. The LSDA density of states in this energy range is 
determined mainly by the covalent Ni-Mn-d hybridiza- 
tion, and by the large exchange splitting of Mn-d elec- 
tron&ii^i^l. At the Fermi level and above a reduced den- 
sity of states is present. The density of states obtained 
from the LDA-I-VGA calculation shows a very strong 
spectral redistribution for the majority spin electrons: 
the LSDA peak situated around — 3eV is lowered in en- 
ergy while in the energy range between -2eV and Ep, 
a spectral-weight transfer towards the Fermi level takes 
place. In particular, the large LSDA-peak at -1.5eV is 
shifted to about — leV, which results in a significant con- 
tribution to the states at the Fermi level. Just above the 
Fermi level, at energies where NQP states are formed in 
the minority-spin channel, a resonance peak is present 
for the majority-spin electrons. A further maximum of 
the density of states is present at 0.5eV. The meaning of 
this maximum will become clear in the Sec. lIII 5] were the 
k-resolved spectral functions are discussed. In contrast 
to our VGA calculation, DMFT results^^ do not change 
significantly the picture for the majority-spin states. Al- 
though the LDA+DMFT density of states shows a similar 
reduction of spectral weight for the peak at — 2eV, its po- 
sition remains unchanged. The differences between these 
two results might be explained by the fact that within 
DMFT Mn and Ni atoms are only coupled via the gen- 
eral many-body and charge-self consistency conditions, 
while correlations are treated independently in the two 
atoms. In contrast, the present VGA approach exactly 



includes correlations on the length scale of the cluster. 
These interatom correlations are possibly responsible for 
the splitting of the covalent Ni-Mn-d electron hybridiza- 
tion in the majority spin states. Due to the breaking 
of this hybridization, the Mn-d exchange splitting is de- 
creased, which could explain the slight shifts of the mi- 
nority unoccupied and occupied majority Mn-d states. 



B. Spectral properties 

In order to gain insight into the nonlocal features of 
the density of states, we compute the k-resolved spectral 
function A(k, w). Majority- and minority-spin spectral 
functions are presented in Figs. [3] and [H respectively, 
with k following high-symmetry points in the Brillouin 
zone (BZ). The explanation of the main features of the 
LSDA band structure was provided by de Groot et. al. 
in his pioneering paper—. Emphasis was placed on the 
interaction between Mn and Sb connected by the sym- 
metry constraint, while less attention was given to the 
Ni atom, although Mn and Ni are first-neighbors and a 
strong hybridization between them is evidenced in the 
density of states. In our LDA-I-VGA calculation, Ni-d 
and Mn-d states are included explicitly, while Sb- states 
are admixed by the downfolding procedure. 




FIG. 3: (Color online) Majority spin LDA-fVCA spectral- 
function of NiMnSb (black/white density plot) along the 
conventional path in the BZ. W(0.5, 1,0), to L(0.5, 0.5, 0.5) 
through r(0, 0, 0), X(0, 1, 0), K{0, 0.75, 0.75) points and end- 
ing at r(0,0,0). The LSDA bands (red, thin solid lines) are 
shown for comparison. Parameters are as in Fig. (2] 

Due to correlation, the majority-spin bands crossing 
the Fermi energy are substantially narrowed with respect 
to the uncorrelated LSDA bands. Specifically, our re- 
sults show for the bands crossing the Fermi level a re- 
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duction of the bandwidth from 3.2eV to 2.2eV. Along 
the path W ^ L both LSDA bands and the VGA spec- 
tral function cross the Fermi level at almost the same 
k-point. The degenerate unoccupied level situated in 
the L-symmetry point, around 1.5eV, is strongly pushed 
towards the Fermi energy, and determines the appear- 
ance of the peak visible at 0.5eV in the DOS discussed 
in Sec. IIII Al At the same time, correlation effects fur- 
ther split the degenerate levels at the F-point seen in 
LSDA at around 2eV. Note that along F — > AT cross- 
ing of the Fermi level occurs close to the corresponding 
crossings of the LSDA-bands. Furthermore, along the 
path X ^ W ^ K both VGA and LSDA bands are 
only weakly dispersive. However, the VGA bands are 
shifted towards the Fermi level, while along the line back 
into the F-point, the Fermi-energy crossing of the VGA 
bands takes place closer to the X-point. 
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FIG. 4: (Color online) Minority-spin LDA+VCA spectral 
function of NiMnSb (black/white density plot) along the same 
BZ path as in Fig. [3] The LSDA bands (red, thin solid lines) 
are shown for comparison. Parameters are as in Fig. [2] 



The minority-spin band structure of LSDA shows an 
indirect gap of about 0.5eV between F and A-point. 
Within this indirect gap formed by the mostly Ai-d occu- 
pied and mostly Mn-rf unoccupied states, the LDA-I-VGA 
results show substantial spectral weight, as can be seen 
in Fig. m Notably, across the Fermi level a weakly- 
dispersive band is present, centered around O.leV, rep- 
resenting the NQP states. At higher energies, in the 
range of 1 to 2eV above Ep, the VGA bands are sub- 
stantially correlation-narrowed with respect to LSDA. 
The features above the Fermi level, including the non- 
quasiparticle states, have dominant Mn-d character. Be- 
low the Fermi level, correlations split off the occupied 
bands having mainly Ni-d character. The spectral weight 
is redistributed: a part is transfered towards the Fermi 



level, however with smaller weight, while most weight is 
transfered towards higher binding energies. The same ef- 
fect is visible in the density of states plot displayed in 
Fig. [51 Notice that while the shift towards higher bind- 
ing energies is also seen in the previous LSDA+DMFT 
calculationi^, the weak shift towards the Fermi level is 
only obtained within the present calculation. 

In order to explore correlation effects in more detail, we 
plot in Fig. [5] the self-energy on Mn sites, in the energy 
range ±2eV around the Fermi level and near the Fermi 
crossing at fc = (0.5, 0.7, 0.3)7r/a. The upper/lower panel 
of Fig. [5] shows the spin resolved real/imaginary part of 
the electronic Mn self-energy. 




E-Ep (eV) 

FIG. 5: (color onhne) Spin resolved self-energy on Mn-d or- 
bitals, at the crossing point k = (0.5,0.7,0.3). 



Below Ep, the minority-spin self-energy is similar to 
the self-energy of the majority spin channel. Just above 
the Fermi level, however, a clear peak in Im(I]yp^) is 
present with a maximum around the energies of the non- 
quasiparticle states (Fig. [5] and Fig. d]). In previous 
DMFT calculations^, a very similar behaviour of the 
imaginary part of the local self-energy was seen. In that 
case, the pronounced feature above Ep was attributed to 
the minority Mn-d{t2g) states. The real part of the self- 
energy displays a negative slope dYi jduj < at the Fermi 
energy for both spin directions, which confirms that the 
quasiparticle weight Z = (1 — is reduced by cor- 

relations. However, while for majority spins |9E|/9a;| is 
clearly less than unity, for minority spins |9E|/9ti;| > 1 
(within our approximation, we cannot determine E with 
sufficient accuracy), suggesting the nonquasiparticle na- 
ture of the minority spin states within the gap. 
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Low-energy spin polarization and comparison 
with experiments 



T" 



• — • Photoemission 
♦— ► P(E) - VCA 
o ^ P(E) - LSDA 




-1.5 

E-E, 



.(eV) 



-0.5 



FIG. 6: (Color online) Energy-dependent polarization ob- 
tained from LSDA and LDA-I-VCA in comparison with data 
from spin-polarized photoemission^. 

To investigate the consequences of the modification of 
majority and minority spectral weight at the Fermi en- 
ergy produced by correlations, we turn to the issue of 
the spin polarization. This is given by the expression 
P{E) = {N^{E) - N^iE))/{N^{E) + N^{E)), N^{E) be- 
ing the spin-resolved density of states, and is plotted in 
Fig. [S] as a function of energy measured from the Fermi 
level. The computed LSDA and VCA values are com- 
pared with the raw data obtained from spin-resolved pho- 
toemission measurement by Zhou et al.^. For this com- 
parison, the density of states was multiplied with the 
Fermi function and a Gaussian broadening of lOOmeV 
was used to account for experimental resolution. Zhu 
et.a^ discuss the appearance of a shoulder close to the 
Fermi level when proper annealing is performed to re- 
store the stoichiometry in NiMnSb. This shoulder is vis- 
ible in the majority spin channel (Fig. 2c from Ref. 
and could be an indication for the correlation-induced 
spectral weight transfer of the majority spin states, not 
present in the LSDA calculations. In addition, the value 



of the spin polarization at the Fermi level obtained from 
our LDA-f VCA-results is situated in the interval of val- 
ues reported exp er iment all y ^ ' ^ '^^ . 

IV. SUMMARY 



We have investigated the effects of correlations in 
NiMnSb using a combined LDA-I-VCA approach. The 
parameters for the effective non-interacting Hamiltonian 
were obtained using the downfolding procedure, for a ba- 
sis including Ni and Mn-d orbitals. The multi-orbital 
Hubbard-type many-body Hamiltonian was solved us- 
ing the Variational Cluster Approach for different val- 
ues of UMn/Ni in the range of 2 — 3eV and JMn/Ni = 
0.65/0.78eV. The results presented do not show signif- 
icant differences for the studied range of parameters, 
nor for different double-counting procedures used. We 
showed that the presence of Ni-d orbitals in the NMTO- 
basis allows for a more complete description of the low- 
energy behavior of NiMnSb. In particular, it correctly 
describes the spectral weight transfer towards the Fermi 
level in the majority spin channel and the formation of 
minority-spin states with vanishing quasiparticle weight 
(NQP states) just above the Fermi level. The analy- 
sis of the minority-spin spectral function shows for the 
NQP states a weakly dispersive band having dominantly 
Mn-d character. Due to electron correlations, the cova- 
lent Ni-Mn d-hybridization in the majority-spin channel 
splits up and part of the weight is transferred towards 
the Fermi level. The simultaneous presence of majority 
spin spectral weight transfer towards the Fermi level, and 
the occurrence of minority-spin non-quasiparticle states 
emphasizes the importance of correlation effects in this 
material, despite the small value of U. 

Despite the fact that high-quality films of NiMnSb 
have been grown, they do not reproduce the half-metallic 
character of the bulk detected by spin-polarized positron- 
annihilation^. On the other hand, one should mention 
that the positron annihilation technique only provides 
an evidence for half metallicity by means of a consistency 
check. In other words, the "proof" is carried out by mod- 
eling the data assuming a half-metallic band structure, 
with a full minority spin gap, from the outsefi''^ . For 
this reason, it would be interesting to revisit the analysis 
of the positron-annihilation data by using the correlated 
band structure obtained here, i. e., by taking into ac- 
count the existence of NQP states. 
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